Dynamic spin susceptibility of superconducting cuprates: 
A microscopic theory of the magnetic resonance mode 
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A microscopic theory of the dynamic spin susceptibility (DSS) in the superconducting state within 
the t-J model is presented. It is based on an exact representation for the DSS obtained by applying 
the Mori-type projection technique for the relaxation function in terms of Hubbard operators. The 
static spin susceptibility is evaluated by a sum-rule-conserving generalized mean-field approximation, 
while the self-energy is calculated in the mode-coupling approximation. The spectrum of spin 
excitations is studied in a homogeneous phase of the underdoped and optimally doped regions. 
The DSS reveals a resonance mode (RM) at the antiferromagnetic wave vector Q = n{l, 1) at low 
temperatures due to a strong suppression of the damping of spin excitations. This is explained by 
an involvement of spin excitations in the decay process besides the particle-hole continuum usually 
considered in random-phase-type approximations. The spin gap in the spin-excitation spectrum at 
Q plays a dominant role in limiting the decay in comparison with the superconducting gap which 
results in the observation of the RM even above Tc in the underdoped region. A good agreement 
with inelastic neutron-scattering experiments on the RM in YBCO compounds is found. 

PACS numbers: 74.72.-h, 75.10.-b, 75.40.Gb 



I. INTRODUCTION 

■ In the superconducting state the spin-excitation spec- 
jtrum of high-Tc cuprates is dominated by a sharp mag- 
netic peak at the planar antiferromagnetic (AF) wave 
Vector Q = 7r(l, 1) which is called the resonance mode 
(RM). It was discovered in the inelastic neutron scatter- 
ing (INS) experiments which revealed a suppression of 
the spectral weight of low-energy spin excitations at low 
temperatures and its transfer to higher energies resulting 
in the RM. There is a vast literature devoted to exper- 
imental and theoretical investigations of spin-excitation 

)rates (a list of references can 
and in a number of publi- 
cations, as e.g., Refs. [6-9]). Here we consider only the 
Plain results of these studies relevant for our purposes. 

[ The RM was discovered at first in the optimally 
doped YBa2Cu30„ (YBCO^) crystal at the energy 
E,- w 41 meV [l3l but later on, the RM was found 
in the underdoped YBCO crystals, Bi2Sr2CaCu208+5 
(Bi-2212) compounds, and other cuprates as well. \2] 
in particular, the RM was observed in the single- layer 



spectra and the RM in cuprai 
be found in the reviews 



cuprate superconductors Tl2Ba2Cu06+2: (Tl-2201) [11 
and HgBa2Cu04+5 (IIg-1201) [l3|, and in the electron- 
doped Pro.88LaCeo.i2Cu04_5 superconductor [l^. This 
demonstrates that the RM is a generic feature of the 
cuprate superconductors and can be related to spin ex- 
citations in a single Cu02 layer. Since the energy of the 
RM was found to scale with the superconducting temper- 
ature. El- w S.S/cbTc in YBCO and Bi-2212 compounds 
and E, w 6 (6.8) fceTc in Tl-2201 (Hg-1201) systems, it 
has been argued that it might constitute the bosonic ex- 
citation mediating superconducting pairing in cuprates 
which has motivated an extensive study of the RM phe- 



nomenon (see, e.g., Refs. 0, [1]). 

The spin-excitation dispersion close to the RM ex- 
hibits a peculiar "hour-glass" -like shape with upward 
and downward dispersions. Whereas the RM energy E^ 
changes with doping, no essential temperature depen- 
dence of £'r and of the upward branch of the dispersion 
has been found. In the optimal doping region the RM 
and both dispersion branches are smeared out above Tc. 
In the strongly underdoped YBCO crystal only the down- 
ward branch is suppressed above Tc, whereas the upward 
dispersion and the RM are observed in the normal pseu- 
dogap state. In particular, a well-defined resonance peak 
at Ey « 33 meV was found in the YBCOe.s crystal in 
the oxygen ordered ortho-II phase with Tc = 59 K at 
hole doping p = 0.09 (Refs. O, Ell). At low tempera- 
ture, T ^ 8 K, the RM revealed a much higher intensity 
than in optimally doped crystals, and it was also seen 
with less intensity even at T 2± 1.4 Tc. 

The extensive study of the twin-free YBCOe.e crystal 
with Tc = 61 K at hole doping p = 0.12 (Refs. [H [T^) 
revealed a different origin of the upward and downward 
parts of the spin excitations close to the RM energy 
Wr = 38.5 meV. Specifically, the high-energy excitations 
do not show noticeable changes at the superconducting 
transition and have the symmetry of the square Cu02 
lattice common to all cuprates. Contrary to this, the low- 
energy part of the spectrum is qualitatively different in 
the superconducting and pseudogap states. In Ref. [l3| 
the INS data was fitted by a spectral function of spin 
excitations in absolute units which was used to calcu- 
late the effective spin-fluctuation-mediated pairing inter- 
action and to estimate the superconducting temperature 
Tc ~ 170 K. [3 

An hour-glass structure of the spin-excitation spec- 
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trum similar to YBCO was found in La2-a;Sr2;Cu04 
(LSCO-a:;) crystals (see Refs. [T9l - [2H ). However, no sharp 
RM was found at T < Tc which may be caused by strong 
disorder effects produced by Sr doping close to the Cu02 
plane or by dynamically fluctuating stripes. A double- 
peaked structure was observed in the local susceptibility 
of the optimally doped LSCO-0.16 crystal ^ with 
Tc = 38.5 K. The lower-energy incommensurate (IC) part 
of the spectrum was found below a; ~ 20 meV which 
may be explained by electron-hole collective excitations 
(or fluctuating stripes), whereas the high-energy upward 
dispersion, peaked at w = 40 — 50 meV, can be viewed 
as overdamped spin-wave-like excitations caused by the 
residual AF exchange interaction J ^ 80 meV. A simi- 
lar, though less pronounced, double-peaked structure in 
x"(aj) was found in the underdoped LSCO-0.085 crystal 
with Tc = 22 K in the pseudogap phase. [2l| 

Generally, the INS experiments suggest the univer- 
sality of the hour-glass structure of spin excitations in 
cuprates, though with material dependent details. The 
high-energy excitations above the RM can be considered 
as overdamped AF spin waves with the exchange inter- 
action decreasing with doping. Since the characteristic 
exchange energy, especially at low doping, is quite large, 
only a weak temperature dependence of the spectrum 
is observed resulting from the increase of the damping 
with temperature that shifts the maximum of the dy- 
namic spin susceptibility (DSS) to a lower energy. The 
low-energy spin dynamics, which reveals a strong tem- 
perature dependence below Tc, may be related to fluctu- 
ating stripe phases with a quasi-one-dimensional order of 
spins and charges or to a nematic (liquid crystal) state 
as discussed later. 

To explain the RM in superconducting cuprates, vari- 
ous theoretical models have been proposed. In a large 
number of studies the Fermi-liquid model of itinerant 
electrons was assumed and the DSS was calculated by 
the random phase approximation (RPA) for the one- 
band Hubbard model or by taking the RPA suscepti- 
bility with the AF superexchange interaction J(q) (see, 
e.g., Refs. [22l - [2^ ) In this approach, the RM is consid- 
ered as a particle-hole bound state, usually referred to 
as a spin-1 exciton. The state is formed below the con- 
tinuum of particle-hole excitations which is gapped at 
a threshold energy Wc < 2A(q*) determined by the su- 
perconducting d-wave gap 2A(q*) at a particular wave 
vector q* on the Fermi surface (FS). In the Fermi-liquid 
approach, the d-wave symmetry of the gap and the shape 
of the FS which should cross the AF Brillouin zone (BZ) 
are essential in explaining the q- and w-dependence of 
the DSS and the RM. 

More complicated electronic models were also consid- 
ered, e.g., in Refs. [6,, .25n28i] , where an extensive study 
of magnetic interactions was performed for the extended 
three-band Hubbard model with large Cu on-site corre- 
lation U. Using the slave-boson representation and 1/N 
expansion (N is the spin degeneracy) the spectrum of 
spin excitations was calculated both for the normal state 



and for the superconducting d-wave state. It was stressed 
that a difference of the magnetic neutron scattering in the 
LSCO and YBCO compounds could be explained by fine 
details of the band structure. (26l . [28| To explain the pe- 
culiarity of the DSS and the RM in the two-layer YBCO 
compounds, a bilayer exchange interaction was invoked 
(see, e.g., Refs. [13, [11]). Usually the downward IC dis- 
persion is well reproduced by fitting the electron interac- 
tion while the intensity of the upward dispersion appears 
to be too weak. To describe the underdoped regime close 
to the insulating (and AF) state, where a model of itin- 
erant electrons cannot be justified, a phenomenological 
spin-fermion model was used (see, e.g., Ref. HI). 

The INS study of the slightly overdoped Bi-2212 crys- 
tal m and of Yi_^Ca^B2Cu306+:E (Y-CaBCO) com- 
pounds [3l| seems to support the spin-exciton scenario. 
In particular, the global momentum shape of the mea- 
sured magnetic excitations is quantitatively described 
within the spin-exciton model with parameters inferred 
from angle-resolved photo-emission experiments on Bi- 
2212 or electronic Raman scattering experiments on Y- 
CaBCO. However, the temperature dependence of the 
RM energy was not studied which should be observed in 
the spin-exciton model due to the temperature depen- 
dence of the superconducting gap. 

The strong correlation limit in the underdoped region 
in cuprates is often treated within the t-J model suit- 
able for consideration of low-energy spin dynamics. To 
take into account the projected character of electron op- 
erators, the slave-boson technique was used (for a re- 
view see Ref. [S^l). In particular, in Ref. [s^ the mag- 
netic excitation spectrum was studied within the t-t'-J 
model using the mean-field slave-particle theory. The 
DSS was calculated in the RPA for spinous in the super- 
conducting and spin-gap states. The RM was revealed as 
the spin-1 exciton below the threshold energy, as in the 
Fermi-liquid models. A qualitative agreement with INS 
experiments on YBCO and Bi-2212 compounds was ob- 
tained. In Ref. 0, a comprehensive analysis of the DSS 
within the mean-field slave-boson theory for the bilayer 
t-J model was performed. A special emphasis was put on 
the explanation of an anisotropy of the spectrum in the 
square Cu02 plane which was related to a d-wave shaped 
FS deformation (dFSD). In the case of a strong dFSD a 
spontaneous breaking of the orientational symmetry of 
the FS can occur leading to its orthorhombic deforma- 
tion (Pomeranchuk instability) . The DSS was calculated 
in the RPA with a renormalized in-plane AF exchange 
interaction. Within the theory, both the collective RM 
and its downward dispersion as well the high-energy ex- 
citations with the upward dispersion were obtained. It 
was also possible to explain the in-plane anisotropy of the 
DSS observed in INS experiments by the dFSD effect. 

However, in the slave-particle theory one has to intro- 
duce a local constraint at each lattice site to reduce the 
enlarged Hilbert space (four states per site) to the phys- 
ical one of the projected electronic states (three states 
per site in the t~J model). In the mean-field approxima- 
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tion the local constraint is relaxed and substituted by an 
averaged one which violates the constraint and leads to 
uncontrollable results. 

A rigorous approach is based on the Hubbard opera- 
tors (HQs) acting in the correct physical space [3j]. In 
Ref. [33, a special diagram technique for the HOs was 
developed and used to calculate the DSS. In particular, a 
generalized RPA was elaborated by summing up different 
types of bubble diagrams. However, the RM excitation 
was not considered there. Calculations of the DSS within 
the HO technique in the conventional RPA in Ref. [32] 
revealed a sharp RM, caused by the opening of the su- 
perconducting gap below Tc, and low-energy collective 
excitations similar to the Fermi-liquid models. 

Let us make a general remark concerning the spin-1 ex- 
citon scenario based on a summation of fermion-bubble- 
type diagrams, as in the conventional RPA. In this sce- 
nario, a strong temperature dependence of the RM en- 
ergy is expected below Tc due to the temperature depen- 
dence of the superconducting gap 2A(q*) and, hence, of 
the threshold energy uJc and E,- < Uc- However, this 
contradicts to experiments. In theoretical calculations, 
usually only the low-temperature limit, T ^ Tc, and the 
T > Tc region are analyzed. It would be important to 
study a region at T < Tc to test the spin-1 exciton sce- 
nario. 

To go beyond a simple RPA in the strong correlation 
limit, the Mori projection technique ^STJ in the equation 
of motion method for the relaxation function has been 
used by several groups (see, e.g., Refs. fssl - lisj ). In the 
t-J model, this method affords to consider the magnetic 
excitations of localized spins in the undoped case within 
the AF Heisenberg model and a crossover to the itin- 
erant electron spin excitations in the overdoped region. 
These studies and, in particular, the DSS calculation and 
the RM analysis in the superconducting state will be dis- 
cussed later and compared with our results. In Refs. [4^ 1 
and [45| . we have formulated a rigorous theory of the 
DSS in the normal state within the projection opera- 
tor method for the relaxation function in terms of the 
HOs. The obtained results, both for the static properties 
(like the staggered magnetization at T = 0, the uniform 
static susceptibility, and the AF correlation length) and 
for the DSS (e.g., the (u;/T)-scaling behavior of the local 
DSS) have shown a good agreement with available cluster 
calculations and neutron-scattering data. In particular, 
it was shown that the kinematical interaction resulting 
from the HO commutation relations plays an essential 
role and gives the major contribution to the damping of 
spin excitations induced by the hopping term t. 

In this respect, we mention a calculation of the DSS us- 
ing the Mori projection technique for the t~J model. [46| 
In this study, the projected character of electron oper- 
ators was neglected by replacing them with the conven- 
tional Fermi operators. As a result, the kinematical in- 
teraction induced by a large hopping term t ~ 3J was 
disregarded and only the exchange interaction J was con- 
sidered. The resulting Mori memory function was repre- 



sented as a sum of two contributions, one determined by 
electrons in the conventional RPA and the second given 
by the spin-spin scattering oc J^. Though the upward 
and downward dispersions of the spin-excitation spec- 
trum and the RM below Tc were reproduced, ignoring 
the hopping term contribution to the memory function 
can give only a qualitative description. 

The discovery of the static IC charge- and spin-density 
waves (CDW and SDW), referred to as stripes, in the 
metallic phase of La2-xBaa;Cu04 (LBCO-x) and LSCO 
compounds has attracted much attention in the recent 
years and was used in the explanation of the IC spin- 
excitation spectrum in hole-doped cuprates (for reviews 
see Refs. [47l-l49j). In particular, in the stripe phase of 
the LBCO-0.125 compound the INS study at T = 12 K 
(> Tc) revealed high-energy spin excitations at w > 
40 meV and a low-energy branch similar to spin waves in 
the two-leg-ladder spin model, Various stripe models 
have been proposed to describe the experimental data. 
In particular, spin-only models for the bond-centered 
stripe with long-range magnetic order in Ref. fsij and a 
two-dimensional (2D) model of coupled two-leg spin lad- 
ders in Ref. [52j were proposed which fitted quite well 
the experimental data. In Ref. |53l |. both bond- and 
site-centered stripes were considered within the time- 
dependent Gutzwiller approximation for the Hubbard 
model. Static stripe spin and charge order coexisting 
with the d-wave superconductivity was studied within an 
extended Hubbard model in the mean-field approxima- 
tion and RPA for the DSS in Ref. [sj]. Using quantum 
Monte Carlo simulations a detailed study of magnetic 
excitations was performed for coupled spin ladders. 55 1 



Magnetic ordering both perpendicular and parallel to the 
stripe direction was found in an array of antiferromagnet- 
ically coupled doped and undoped two-leg ladders. (5a | 

However, the static stripes have not been detected 
in moderately doped YBCO crystals, where only dy- 
namically fluctua ting stripes can be suggested (see, e.g., 
Refs. d, [T3 - [l6l [571). Moreover, instead of a rigid 
stripe array proposed in Ref. [HI], the 2D character of 
the IC spin fluctuations was observed in the untwinned 
YBCOe+x crystals with x = 0.6 and 0.85. [s^l The IC 
peaks exhibiting the in-plane anisotropy were explained 
as a result of a possible liquid-crystalline stripe phase or 
a nematic state. Therefore, though the inhomogeneous 
stripe picture can be applied for LSCO and LBCO sys- 
tems, it seems to be not a universal origin of the hour- 
glass spin-excitation dispersion and the RM in cuprate 
superconductors. Inhomogeneous phases of cuprates re- 
lated to fluctuating stripes with a quasi-one-dimensional 
order of spins and charges or to a liquid crystal state may 
be important only in the explanation of the low-energy 
collective spin fluctuations. 

To clarify some of the open problems in describing 
the RM phenomenon, such as the appearance of the 
RM above Tc and its weak temperature dependence, 
in the present pager we extend our microscopic theory 
(Refs. 44g and [45[) to the superconducting state. Al- 
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though our general formulation for the DSS is similar 
to the original Mori memory function approach used in 
other publications, as in Ref. [s^, in the previous studies 
of the t-J model only the bubble-type diagrams simi- 
lar to the RPA were considered which ignores the im- 
portant role of spin excitations in the decay process. 
The energy gap at the AF wave vector Q of the or- 
der of the RM energy in the spin-excitation spec- 
trum strongly reduces the damping at low temperatures, 
T <^ ^ b k^Tc, which results in the emergence of a 
sharp peak in the spectral function. In the low doping 
region, where the damping is extremely small, the RM 
is found even above T^. In the overdoped region, at hole 
concentration S ^ 0.2 and high Tc, the spin-excitation 
damping becomes large and the opening of the supercon- 
ducting gap enhances the intensity of the RM, so that 
it becomes observable only below Tc. So, as compared 
with the spin-exciton scenario, we propose an alternative 
explanation of the RM and the upper branch of the dis- 
persion which are driven by the spin gap at Q instead 
of the superconducting gap 2A. A good agreement of 
our results for the temperature and doping dependence 
of the spin-excitation spectrum and the RM with INS 
experiments provides a strong support for the proposed 
theory. We have not found a lower branch of the hour- 
glass spectrum which may be related to inhomogeneous 
states in the Cu02 plane neglected in our analysis. 

In the next section we present the basic formulas for 
the DSS and the self-energy which are a generalization of 
our theory in Ref. [45|] to the superconducting state. Nu- 
merical results for the spin-excitation spectra are given in 
Sec. Ill, where the temperature and doping dependence 
of the damping and the RM are discussed. The conclu- 
sion is given in Sec. IV. Details of the calculations within 
the mode coupling approximation (MCA) are presented 
in the Appendix. 



II. RELAXATION-FUNCTION THEORY 
A. Dynamic spin susceptibility 

It is convenient to consider the t-J model in the Hub- 
bard operator representation 



(7(7 \ 
3 ) 



(1) 



interaction. The Hubbard operators A^ 



describe transitions between three possible states at a 
site i on a square lattice: an empty state \i,a) — |i, 0) 
and a singly occupied state \i,a) = |i,(T) with spin a = 
±(1/2), (ct = — (t). The number and spin operators in 



terms of the Hubbard operators read: 



517 



(2) 



The Hubbard operators obey the completeness relation 
+ J2a ~ 1 which preserves rigorously, con- 

trary to the slave-boson approach, the constraint of no 
double-occupancy of any lattice site. The Hubbard op- 
erators have the commutation relations 



7/3 



which results in the kinematical 



interaction. Here, the upper sign pertains to Fermi-type 
operators like X^" changing the number of electrons, and 
the lower sign pertains to Bose-type operators, such as 
the number operator or the spin operators, Eq. The 
chemical potential /i is determined from the equation for 
the average electron density n = {Ni) = 1 — 5, where 
5 = (A™) is the hole concentration. 

In Ref. [13], applying the Mori- type projection tech- 
37[ , elaborated for the relaxation function, we have 



nique 



derived an exact representation for the DSS x{^t'^) re- 
lated to the retarded commutator Green function (GF) 
(see Ref. Hi), 



x(q,c.) = -((5+|5Zq)). = 



m(q) 



- uj S(q, cj) — a;2 ' 



(3) 



where TO(q) = {[iS+, SZ^]) - ([ [5*+, if], SlJ), and c^q is 
the spin-excitation spectrum in a generalized mean-field 
approximation (GMFA). The self-energy is given by the 
many-particle Kubo-Mori relaxation function 



S(q,c.) = [l/m(q)]((-5+|-^ 



(propcrj 



(4) 



where -S^ = [[S^,H], H] (for details see Ref. ;44i). 
The Kubo-Mori relaxation function and the scalar prod- 
uct are defined as (see, e.g., Ref. [60[ ) 



/>oo 

{{A\B))^ = -i dte''^'{A{t),B), 
Jo 



(5) 



and 



iA{t),B) 



dX{A{t ~iX)B), (i^l/k^T, (6) 



respectively. The "proper" part of the relaxation func- 
tion dl]) does not contain parts connected by a single 
zero-order relaxation function which corresponds to the 
projected time evolution in the original Mori projection 
technique [S^. The spin-excitation spectrum is given by 
the spectral function defined by the imaginary part of 
the DSS (121), 



where is the hopping integral and Jij is the exchange x'{^.^'^) 



-u S"(q, io) m(q) 



-L^I]'(q,c^)]2 + [c^S"(q,w)] = 



(7) 



where i;(q,w + iO+) = T.'{q,u}) + il]"(q,cj), and 
T.'{q,uj) = -S'(q,-w) and S"(q,w) ^ S"(q, -w) < 
are the real and imaginary parts of the self-energy, re- 
spectively. 
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B. Static susceptibility 



C. Self-energy 



The general representation of the DSS ([3]) determines 
the static susceptibiUty Xq — xili 0) by the equation 

Xq=(5+,5Zq) = m(q)/c.2. (8) 

To calculate the spin-excitation spectrum ujq the equality 

m(q) = (-5+, 5Zq) = iol {S+, SZq), (9) 

is used, where the correlation function (— S'+,5Zq) is 
evaluated in the GMFA by a decoupling procedure in 
the site representation as described in Ref. [i^. This 
procedure is equivalent to the MCA for the two-time 
correlation functions. This results in the spin-excitation 
spectrum 



ujI = 8t2Ai(l-7q)(l-n-i;^2,o-2Fi,i) 
+ 4j2(l-7q)[A2| -aiCi^o(47q + l) 
+ a2(2Ci.i + C2,o)], (10) 



where t and J are the hopping integral and the exchange 
interaction for the nearest neighbors, respectively, and 
7q = (1/2) (cos -I- cos Qy) (we take the lattice spacing 
a to be unity). The static electron and spin correlation 
functions are defined as 



R = (^o°^r) = ^E^i<5"'''' (11) 
q 

lECqe^'^^ (12) 
q 



where R = nax + may . The GMFA spectrum pO|) is 
calculated self-consistently by using the GMFA approxi- 
mation for the static correlation function p^ . 



(13) 



Cq=!!^COth^ 



The decoupling parameters ai, a2 and Ai, A2 in Eq. (|TOl) 
take into account the vertex renormalization for the spin- 
spin and electron-spin interaction, respectively, as ex- 
plained in Ref. (45| . In particular, the parameters a 1,02 
are evaluated from the results for the Heisenberg model 
at (5 = and are kept fixed for S ^ 0. The parameters 
Ai, A2 are calculated from the sum rule Co,o = {Sq Sq) — 
(1/2)(1 - S) with a fixed ratio A1/A2 =' 0.378. In the 
superconducting state, the electron correlation function 
Fr is calculated by the spectral function for electrons in 
the superconducting state (see Eq. (|22|) V The variation 
of the parameters Ai, A2 below the superconducting tran- 
sition is negligibly small and practically has no influence 
on the spectrum ojq. 

The direct calculation of m{q) yields 



m(q) = -8 (1 - 7q) [tFi^o + ^t^i.o] 



(14) 



Thus, the static susceptibility (jS]) is explicitly determined 
by Eqs. pll |) and pl |) . 



In what follows, we consider the t-J model at a finite 
hole doping S > 0.05 when, as discussed in Ref. (45| , 
the largest contribution to the self-energy is Et(q, w) 
coming from the spin-electron scattering. It is deter- 
mined by the hopping term Ht in the t-J model ac- 
cording to the equation for the spin-density operators: 
— S'q = [ [S^,Ht], Ht]. As described in the Appendix, in 
the MCA this contribution reads 



S;'(q,c.) = - 



7r(2t)4(e'9-' - 1) 



dujiduj2duj3 



m(q) uj 



qi:q2 



Sq.(-2)[(Aq.q..q3+Aq3,q..qJ<(-lX3(-3) 



qi.q2,q3 ' q3,q2,qi'' qi^ 

<,(a.i)A^3,(.;3)], (15) 



2Aqj^ q2.q3 Aqg q^^qj^ ^iq^ 



where q3 = q— qi— q2. The Fermi and Bose functions are 
denoted by n(w) = (e^'^ + and N{lo) = ( 



The vertex function A, 



qi,q2,q3 



= (e/3" - 
is defined by Eq. 



1 



Here we introduced the spectral functions: 



q.(^) 



-(l/7r)Im((X^n^q")) 



q 

05- \ 



-{l/^)lu.{{Xl^\X%)), 
Bq{u) = (l/^)x"(q,c.). 



(16) 
(17) 
(18) 



where A^'^ (uj) are determined by the retarded anticom- 
mutator GFs for electrons (see Ref. [1^). In comparison 
with the expression for the self-energy in the normal state 
considered in Ref [i^, in Eq. (fT5|) there is the contribu- 



tion proportional to the anomalous GF ((Xq'^|X'lq))^ 
which is nonzero in the superconducting state only. 

It should be emphasized that the self-energy (fT5|) is de- 
termined by the decay of a spin excitation with the energy 
Lu and wave vector q into three excitations: a particle- 
hole pair and a spin excitation. This process is con- 
trolled by the energy and momentum conservation laws, 
UJ = (wa — uji) + UJ2 and q = qi + q2 + Qs, respectively. 
In previous studies of the t-J model the contribution of 
the additional spin excitation has been neglected (see, 
e.g., Ref. [11]) or approximated by static or mean-field- 
type expressions (see, e.g., Refs. 38 1 and [i^l)- That 



is, in these approximations the spin-excitation contribu- 
tion was "decoupled" from the particle- hole pair. We 
can derive the particle-hole bubble approximation from 
Eq. (jlSp , if we ignore the spin-energy contribution uj2 in 
comparison with the electron-hole pair energy, or, equiv- 
alently, if in the MCA, Eqs. ([Ml) and (KTl . the time- 
dependent spin correlation function is approximated by 
its static value: {SZqS+{t)} ~ (S'ZqS'+) = Cq. More- 
over, excluding the spin-excitation wave- vector q2 from 
the wave- vector conservation law, we have q = qi + qs. 
As a result of these approximations in Eq. (jl5l) . we ob- 
tain the self-energy in the form of the particle-hole bubble 
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approximation: 



S'/(q,^) 



duji [n{uji) — n{LUi + co)] 



N 



qi 



'^qi,q-qi Kl'^i^^) ^q-qi^('^l + ^)] ' 



(19) 



where the averaged over the spin-excitation wave vector 
q2 vertexes are introduced, 



qiqa 



]^ E [^qi.q2,q3 + 



■p^ / , ^q2^qi,q2,q3^q3,q2,qi- 
q2 



, (20) 
(21) 



In the approximation (jl9p only the opening of a supercon- 
ducting gap in the particle-hole excitation can suppress 
the damping of spin excitations due to the decay into 
particle-hole pairs which may result in the RM. Below we 
discuss in more detail why a particle-hole bubble approx- 
imation for the self-energy, Eq. p^ . leads to a different 
behavior of the spin-excitation damping in comparison 
with the results obtained for the full self-energy ([T5|) . 



III. RESULTS AND DISCUSSION 



A. Self-energy approximation 



In the calculation of the self-energy dlSt we adopt the 
mean-field approximation (MFA) for the electron spec- 
tral functions (jl6l) and (jl7p which in the superconducting 
state can be written as 



^1 + £q 

2^1 



d{co-LOi), (22) 



AU^) = Q E 



A 



CJl=±-Ec 



2uJi 



S{U! — UJl) . 



(23) 



Here Q = 1 — n/2 is the Hubbard weighting factor and 
the superconducting gap function Aq^- = (sgn a) Aq. 
In the electron spectrum Sq we take into account only 
the nearest-neighbor hopping t and consider the en- 
ergy dispersion in the Hubbard-I approximation: £q = 
—At Q 7q — /i. The spectrum of quasiparticles in the su- 
perconducting state is given by the conventional formula 

Eq = ^Sq + A^. For the spin-excitation spectral func- 
tion (fT8|) we take the form: 

i?q(a.)-m(q) ^-^(^-^i)' (24) 

where the spectrum of spin excitations ujq is determined 
by the pole of the DSS, Wq = [w^ Wq S'(q,u;q)]i/2 . 



Here, the real part of the self-energy E'(q, w) is calcu- 
lated perturbationally by taking the spectral function 
(f24| with the GMFA spectrum wi = iwq. Using these 
spectral functions, after integration over the energies in 
Eq. (|15p we write the imaginary part of the self-energy in 
the following form convenient for the calculation in the 
limit T ^ 0: 



S;'(q,c.) 



7r(2t)4 
a;m(q) 



E E E E 



m(q2) 
[(A^ 



qi,q2 toll=±i5qj IjJ2=±I^ci2 "3 = ±-Eq 

N{uj2)n{—uji)n{uj3) + N{-u}2)n{uJi)n{—uj3) 



+ A2 

qi.q2,q3 ' q3,q2,qi 



8UJ1UJ2UJ3 

){UJI + Sq,){LJ3 + Sq,) (25) 



-2Aq, 

,q2,q3Aq3,q2,qi^qi^q3] ^{^ + ^^l 1^2 ~ ^s)- 

Similar calculations for the self-energy in the 

particle-hole bubble approximation yield 



S;'(q,c.) 



7r(2t) 



4 r,2 



■m{q)uj N 



EE E 

4^^^^ K,q-qi(^l + £qi)(^2 + eq-qj 

-^qi.q-qi^qiN-qJ + ^1 " ^2). (26) 

We consider both the d-wave and s-wave symmetry 
of the superconducting gap which we write as Aq'*' — 
(A/2)(cosga; — cosqy) and A*^''-' — A with the tem- 
perature dependent amplitude A(r). In mmierical cal- 
culations we assume that A(T) follows the conven- 
tional Bardeen-Cooper-SchriefFer (BCS) theory. In par- 
ticular, A(T)/fcBTc = 1.76, 1.72, 1.6, 1.24 for T/T^ = 
0, 0.4, 0.6, 0.8, respectively. By taking, instead of the 
BCS ratio 2Ao/kBT^ = 3.52, Aq = A(T = 0), the ra- 
tio 2Ao/A:bT'c = 4.3 for a pure d-wave superconductor 
(see, e.g., Ref. [6l|) we have found that the results do 
not change noticeably, i.e., less then 5% at T = and 
even less at finite temperatures. We mainly consider 
two doping values, 6 — 0.2 which is larger than opti- 
mal doping and 6 — 0.09 for the underdoped case. For 
S = 0.2 we fix the superconducting transition temper- 
ature as fceTc — 0.025t, while for 6 = 0.09 we take 
/cbTc — 0.016t. For the hopping parameter t — 0.313 eV 
these values are close to Tc = 91 K in the nearly optimally 
doped YBC06.92 single crystal in Ref. [J and Tc = 59 K 
in the underdoped {S = 0.09) YBCOe.s crystal studied 
in Ref. [l3|- We take the exchange interaction J = 0.3 1 
and measure all energies in units of t . 



B. Spin-excitation damping 

To elucidate the role of spin excitations in the damp- 
ing and their relevance to the shape of the spectral func- 
tion, Eq. (O, we consider the temperature dependence 
of the spin-excitation damping at the AF wave vector. 
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(o/t 

FIG. 1: Spin-excitation damping r(Q,aj) for S — 0.2 at T < 
Tc for (a) the d-wave and (b) s-wave pairing, and (c) in the 
normal state. 



r(Q,cj) = -(l/2)S;'(Q,w). Figure [U shows the damp- 
ing for 6 = 0.2 at various temperatures in the case of 
the d-wave (a) and s-wave (b) pairing. The difference of 
the damping appears only at low uj and T. In particular, 
the damping for the s-wave gap at T = disappears at 
ui < 2Aq ~ 4fcB'7c = O.lt, while for the d-wave gap it van- 
ishes at w ~ Aq ~ 0.05t. A weak damping was obtained 
also for the normal state shown in Fig. [1] (c) when the 
contribution from the superconducting gap functions in 
the self-energy I]"(Q,a;), Eq. is omitted. The sim- 
ilar smooth variation of the damping with energy in all 
three cases below Tc, contrary to a step-like dependence 
obtained in the particle-hole bubble approximation (see 
below) , demonstrates that the superconducting gap plays 
a minor role in suppressing the damping, and the gap ojq 
in the spin-excitation spectrum is responsible for such a 
peculiar behavior. 

For a lower doping, the damping becomes an order 




0.20 



FIG. 2: Spin-excitation damping r(Q, oj) for 5 = 0.09 for the 
d-wave pairing at T = (solid line), T = Tc (dotted line), and 
T = lATc (dash-dotted hne). The damping at T = for the 
electron dispersion with t' /t — 0.37 and t" /t — 0.1 is shown 
by the dashed line. 



of magnitude weaker, as shown in Fig. [5] at ^ = 0.09, 
even above Tc (T — lATc). For comparison with re- 
sults obtained in the fermion-bubble approximation, we 
calculate the damping also for a more general electron 
dispersion, Sq = —At 7q + At' j'^ — At" 7q — fj,, taking into 
account hopping between next- and third-nearest neigh- 
bors, t' and t" , respectively, where 7q = cosg^ cosgj, and 
7q — (l/2)(cos2g^ -I- cos2qy). In Fig. [5]the damping at 
T = for the parameters t'/t = 0.37 and t" /t = 0.1, 
proposed in Ref. 18| for the electron dispersion in the 
antibonding band in YBCOg.e, is shown by the dashed 
line. The difference between the damping calculated for 
the electron dispersion gq = —At Q 7q— /i in Eqs. ([^ and 
and the more general one with nonzero parameters t' 
and t" is negligible. This proves that the shape of the FS 
and the electron dispersion are unessential factors in our 
theory. To describe the electron-doped cuprates within 
the t-J model, we should change the signs of the hopping 
parameters tij in Eq. ([l}, in particular, —t^t>0. This 
replacement does not change the obtained results and, 
therefore, the theory can be also applied to the electron- 
doped cuprates with a proper fit of model parameters. 



Sometimes, the RM observed above Tc in the under- 
doped cuprates is related to a pseudogap in the elec- 
tronic spectrum (see, e.g., Ref. Q). We propose an- 
other explanation: The RM above Tc is the result of 
the gapped spin excitations in the self-energy (1^ lead- 
ing to a very weak damping in the underdoped region 
which is outlined in more detail below. This explanation 
is supported by studies of the spin-excitation damping 
Fq = -(l/2)S;'(q,cj = Wq) at T = shown in Fig. H 
The small difference between the damping in the d-wave 
superconducting state and the normal state observed for 
the full self-energy, Eq. ([25]), confirms that the supercon- 
ducting gap does not play an essential role in suppressing 
the damping Fq , in particular in the underdoped region. 
At the same time, the sharp increase of Fq away from 
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FIG. 3: Spin-excitation damping Fq for (a) S — 0.2 and for 
(b) S = 0.09 at T = for the d-wave pairing (solid line) and 
in the normal state (dashed line). 
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FIG. 4: Spectral function x"(Qi for the d-wave and s-wave 
pairing in comparison with the normal state at T = OATc for 
S = 0.2. 



the AF wave-vector Q explains the resonance character 
of spin excitations at Q. 

Though the damping in Fig. [1] looks similar, the spec- 
tral functions shown in Fig. |4] for (5 = 0.2 at T = 0.41; 
reveal a strong enhancement of the intensity of the RM in 
the superconducting state. In comparison to the normal 
state, where the contribution from the superconducting 
gap is omitted, the peak intensity is about two (five) 
times larger for the d (s)-wave symmetry of the gap. 

Quite a different behavior of the damping and the 
spectral function is obtained for the reduced self-energy, 



Eq. (|26p . with a contribution only from a particle- hole 
bubble. Figure [5] shows our results for the spectral func- 
tion ^"(ctjw) and for the damping F(Q,a;). To com- 
pare these functions with those calculated in Ref. [s^ , we 
adopt the electron dispersion used in Ref. [11], s^-l" = 
4te//7q — At'^j: r cosqx cosqy with t^ff = 0-3t and 



eff 



eff 

-0.lt and take the gap parameter Aq = O.lt. 
The obtained results are quite close to those shown in 
Fig. 1 of Ref. m (where x'L{<l,^) = (l/2)x"(q, c^) is 
plotted). At T = 0, we observe a much narrower RM, 
but with a lower intensity in comparison with the RM 
calculated with the full self-energy, Eq. (|25p . as shown 
in Fig. |6l The energy E'r of the RM shown in Fig. [5] (a) 
noticeably decreases with increasing temperature, con- 
trary to a negligible shift of the RM shown in Fig. [5] for 
T — 0.4rc. This comparison demonstrates that in the 
particle-hole bubble approximation the superconducting 
gap plays a crucial role in the occurrence of the RM with 
E,{T) < 2A(r), while in the full self-energy ^ the su- 
perconducting gap and details of the electron dispersion 
are less important. In particular, for the underdoped case 
S = 0.09 we have not found visible changes of the damp- 
ing function shown in Fig. [5] for the electron dispersion 
with t' = and t' = —O.lt. For the reduced self-energy, 
Eq. (PB)) . the damping vanishes for both types of the dis- 
persion in the underdoped region, and in order to obtain 
a finite damping at 6 = 0.1 in Ref. [s^ (see Fig. 2), the 
authors have to use the electron density of states in the 
damping function (see their Eq. (20)) instead of the q- 
dependent electron spectral functions. 

This difference can be explained as follows. Whereas in 
the particle-hole bubble approximation given by Eq. ([26]) 
the spin excitation with the energy lo at the wave vec- 
tor Q can decay only into a particle-hole pair with the 
energy uj{Q) = -Eq+q + -Eq, in a more general process 
described by Eq. (|25p an additional spin excitation par- 
ticipates in the scattering. In the limit T ^> 0, the de- 
cay process is governed by another energy-conservation 
law, uj{Q) — Eq^ + Eq^ + Wqj, whcre the largest contri- 
bution from the spin excitation comes from ujq^ ~ ujq 
due to the factor rn(q2), Eq. ([HI), in Eq. (gHl)- This 
energy-momentum conservation law strongly reduces the 
phase space for the decay and suppresses the damping 
of the initial spin excitation with the energy w(Q). In 
fact, the occurrence of an additional spin excitation in 
the scattering process with the finite energy wq plays 
a role similar to the superconducting gap in the exci- 
tation of the particle- hole pair in Eq. ([25)1 . Therefore, 
the damping at low temperatures {k^T ujq ^ E^) 
appears to be small even in the normal state as demon- 
strated in Fig. [1] (c). In the case of the particle- hole 
relaxation, the condition for the occurrence of the RM, 
w(Q) = Eq+Q + Eq < 2A(q*), imposes a strong re- 
striction on the shape of the FS which should cross the 
AF Brillouin zone to accommodate the scattering vector 
Q and the vector q* on the FS. In the case of the full 
self-energy, Eq. (j25p . the energy- momentum conservation 
law for three quasipartricles does not impose such strong 
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FIG. 5: (a) Spectral function x"(Qi ^) ^nd (b) spin-excitation 
damping r(Q,ix)) calculated in the particle-hole bubble ap- 
proximation, Eq. (|19|) . at 5 = 0.2 for the d-wave pairing 
(Ao = O.lt taken from Ref. [H ) at T = (solid line) and 
T = 0.4Tc (dotted line), and for the normal state at T = 
(dashed line) and T = Tc (dash-dotted line). 



limitations. 



C. Resonance mode 

Experimentally, the RM energy decreases with un- 
derdoping following the superconducting transition tem- 
perature, Et: ~ 5.3fcBTc, but only weakly depends on 
temperature (see, e.g., Refs. 0, Now we discuss the 
temperature and doping dependence of the RM and its 
dispersion within our theory for the d-wave pairing. 

The temperature dependence of the spectral function 
in the overdoped case 5 = 0.2 is shown in Fig. |6l It 
has high intensity at low temperatures, but strongly de- 
creases with temperature and becomes very broad at 
T '--^ Tc as found in experiments (see Ref. (H). In Fig. [7] 
the temperature dependence of the spectral function for 
the underdoped case 5 = 0.09 is plotted. Whereas the 
resonance energy decreases with underdoping, the inten- 
sity of the RM greatly increases in accordance with exper- 
iments. The RM energy weakly depends on temperature 
and is still quite visible at T = Tc and even at T = ^AT^. 

The dispersion of the spectral function for 5 = 0.2 
is shown in Figs. |8] and [9l A strong suppression of the 
spectral-function intensity away from Q = 7r(l, 1) even at 
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FIG. 6: Temperature dependence of the spectral function 
x"(Q,t^) at 5 = 0.2. 
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FIG. 7: Temperature dependence of the spectral function 
x"(Q,t^) at 5 = 0.09. 



T — Q explains the resonance-type behavior of the func- 
tion at low temperatures. This suppression of the inten- 
sity is in accord with the sharp increase of the damping 
away from Q shown in Fig. [3l As can be seen from 
Figs. [8] and [U we obtain an upward dispersion of the 
resonance energy f2(Q -I- q). Our numerical results may 
be well fitted by the quadratic dispersion law for small 
wave vectors, (qx,(ly) ^ 0.27r , 



n{Ql + Cl) = [El+c\ql+ql)] 



1/2 



(27) 



where Ejt = 0.12(0.09) and c/at ^ 0.36(0.45) for 
5 — 0.2 (0.09). In the conventional units we obtain 
E,, = 38 meV, c = 425 meV-A for 5 = 0.2 and 
Er = 28 meV, c = 545 meV-A for 5 = 0.09, where 
we take a — 3.82 A and t — 313 meV. Note that the 
quadratic dispersion was reported in several papers. For 
example, the dispersion p7)) was found in Ref. [1^ for the 
acoustic mode in YBCOe.s {5 = 0.09) with the parame- 
ters Ey. = 33 meV and c ~ 360 meV-A that qualitatively 
agrees with our results. 

We have not found the downward dispersion at en- 
ergies below the RM detected in neutron-scattering ex- 
periments, for exarnple, on YBCOg.s (Ref. [iH]) and on 
YBCOe.e (Refs. However, as argued in Ref. [H, 
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FIG. 8: Spectral function x"(q,'^) for the wave vectors (a) 
q = 7r(l, ^) and (b) q = 7r(^, ^ at T = for 5 = 0.2. 
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FIG. 10: Real part of the self-energy E'(Q,ctj) for the d-wave 
pairing in comparison with the normal state at T = for 
S = 0.2. 
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FIG. 9: Spectral function x"{^,^) near the wave vector Q : 
7r(l, 1) at T = for (5 = 0.2. 



two distinct regions of spin excitations may be suggested: 
a low-energy part below w 33 meV, which can be de- 
scribed as incommensurate stripe-like collective spin ex- 
citations for the acoustic mode, and a high-energy part, 
which has a spin- wave character. The high-energy part of 
the spin excitations has an isotropic in-plane dispersion 
while the low-energy excitations show a one-dimensional 
character. The different nature of the two parts of the 
spectrum is also revealed in their temperature depen- 
dence: the low-energy acoustic part of the spectrum is 
strongly influenced by the superconducting transition, 
while the high-energy part of the spectrum does not 
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FIG. 11: Energy of the resonance mode i?r at T = as a func- 
tion of hole doping S in comparison with experimental data 
for YBCO from Ref. (open circles), Ref. H (full squares), 
and Refs. [l|, [3 (fuU circles). 



change appreciably with temperature up to 85 K. Sim- 
ilar differences between the upward and downward dis- 
persions were found in Refs. [ll] and [l7j . 

The downward dispersion of the lower-energy part of 
the spectrum was explained within the stripe-like mod- 
els, as discussed in the Introduction. It was also found 
for a Fermi liquid model in the RPA approach (see, e.g., 
Ref. [131) or, for the t-J model, within the particle-hole 
bubble approximation (see Refs. [1^, SO, IH | ) . The dis- 
persion was explained by a special wave-vector depen- 
dence of the particle-hole bubble diagram related to the 
wave- vector dependence of the dx^^y2 superconducting 
gap and to a specific for cuprates 2D FS. Since in our 
theory beyond the RPA the RM energy does not crit- 
ically depend on specific properties of the FS and the 
superconducting gap, the downward dispersion cannot 
be found. To discuss this problem in detail, the role of 
stripe excitations in the spin-excitation spectrum should 
be elucidated. 

Now we discuss the doping dependence of the RM en- 
ergy E^{5). At low temperatures, the real part of the 
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FIG. 12: Spectral function x"(Q,'^) for doping^ = 0.2 com- 
pared to experimental data for YBCO6.92, Ref. [J], at T = 5K 
(squares) and T = lOOK (circles). 
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FIG. 13: Spectral function x"{Q,'^) for doping 5 — 0.09 com- 
pared to experimental data for YBCOe.s, Ref. [14|, at T = 8K 
(squares) and T — 85K (circles). 



self-energy I]'(Q,a;) < is quite large as shown in 
Fig. [TUl This considerably softens the energy of spin ex- 
citations wq in the GMFA, Eq. pUj) . shifting the pole 
of the spectral function x"(Qi'^) to a lower energy: 
WQ = [w^-WQ |E'(Q,wq)|]^/^ . Experimentally, the RM 
energy is measured by the position E,- of the maximum 
in the spectral function x"(Qj'^) which deviates from 
ujQ due to a finite width of excitations. In Fig. [TT] the 
doping dependence E,-{6) in the superconducting state 
at r = determined by the maximum of the spectral 
function is plotted. Thereby, for the doping dependence 
of Ao(5) = 1.76Tc(i5) we used the universal empirical 
formula T,{d) = Te,„,ax [1 - ^ (S - V)'] ^^ere 
Sopt = 0.16, Tcmax — 93 K, and the value of /? = 75 
was fitted to obtain Tc = 59 K for S = 0.09 in YBCOg.s, 
Ref. [13] • With decreasing 5, decreases which quali- 
tatively agrees with the experimental data. The energy 
of the RM tends to zero at the critical doping 5c = 0.038 
below which the long-range AF order emerges at T = 0, 
as we have shown in Ref. |45| . So, in our scenario the RM 
is just the soft mode which brings about the long-range 
AF order below the critical doping. 

Experimentally, in the overdoped region the RM en- 
ergy decreases with increasing doping (see Fig. [TT] and, 
e.g., Ref. [13]), while in our theory E^ tends to increase 
due to the increasing energy ujq. The agreement with 
experiments in Fig. [TT] can be improved in the optimally 
doped region, 5 < 0.2, if we change the parameter t to 
t = 0.4 eV instead of the adopted value t = 0.313 eV. 
We also note that in the approach of Ref. [38| a too large 
superconducting gap Aq ^ O.lt (in comparison with our 
value Aq ^ 0.044 1) has to be taken to fit the RM energy 
to the experimentally observed one. 

In Fig. [12] we compare our results with the neutron- 
scattering data for the nearly optimally doped YBCO6.92 
single crystal [l\ at T ^ 5K and T = lOOK. In this 
sample, Tc = 91 K and the RM energy E^ ~ 40 meV= 
5.1fcBrc > 2Ao (taking 2Ao((5) = 3.52fcBTc(5) we have 
Er ~ 2.9A0). For 6 = 0.2, our calculations yield 
E, = O.m = 38 meV= 4.8fcB Tc = 2.7 Aq (t = 0.313 eV, 



/cB Tc — 0.025t; see Sec. III. A). Concerning the heights 
of the peaks, the scale of the experimental data, given 
in arbitrary units, was adjusted by fitting only the peak 
height at T = 5K to our result at the same temperature. 

In Fig. [13] our results are compared with the experi- 
mental data for the underdoped ortho-II YBCO6.5 single 
crystal with E^ = 33 meV= 6.5 ksT^ = 3.7 Aq at T = 8K 
and T = 85K (see Fig. 14 in Ref. For 5 = 0.09, our 

theory gives E, = 0.09t = 28 meV= 5.6 ksTc = 3.2 Aq. 
Here, the experimental peak heights were also scaled by 
fitting only the peak height at T = 8K to our result at 
the same temperature. We note a weak temperature de- 
pendence of the RM energy observed experimentally and 
obtained in our calculation. In both compounds the RM 
energy is larger than the superconducting excitation en- 
ergy, 2Ao, while in the spin-1 exciton scenario the RM 
energy E,- has to be less than 2Ao. So we obtain a good 
agreement of our theory with neutron-scattering experi- 
ments on YBCO crystals both near the optimal doping 
and in the underdoped region. 

It is also interesting to compare our results with the 
recent detailed experimental study of the spin-excitation 
spectrum reported in Refs. [l^ and [13] for the twin-free 
YBCOe.e crystal with Tc = 61 K at hole doping S = 0.12. 
For that, we calculate the spectral function x"(q, w) for 
the same hole doping and Tc. From the experimental 
value Wr — 38 meV and damping F = 11 meV we cal- 
culate i?r = 36.4 meV [s^. To get agreement with the 
theoretical value E^ = 0.091<: we adopt t = 0.4 eV, as 
also inferred from the ARPES data in YBCOe.e- [H For 
the upper branch of spin excitations above the RM we 
obtain a quadratic dispersion given by Eq. (|27|) . whereas 
m Ref. [17], Eq.(ll), the dispersion was fitted by a quar- 
tic power law. 

For a more detailed comparison of the RM disper- 
sion with the experimental data in Ref. [l3] for the up- 
per branch, in Fig. [14] we present the spectral function 
X"(q,w) at T = and T = 70 K for the wave vec- 
tors q — TT{qx,l) at energies ui/E,: = 1, 1.24, and 1.58 
which, putting E,- = uj,-, corresponds to a; = 38, 47, 
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1.00 



CT 200 



1.00 



FIG. 14: Spectral function x"(ci) ^) for <l ~ 1) at doping 

5 = 0.12 and energies ui/Er — 1 (bold line), 1.24 (dashed line), 
and 1.58 (dotted line) at (a) T = and (b) T = 70 K (in units 
of t = 0.4 eV). 



and 60 mcV in Ref. [17[. The damping of the RM in 



our theory stron gly depends on q as shown in Fig. |31 
whereas in Ref. |17| the damping was taken to be q- 
independent, F = 11 meV. In our calculations we get 
a much smaller damping: Fq = 0.0034t ~ 1.4 meV at 
T = 0. These differences result in much sharper spectral 
functions in our theory in comparison with the experi- 
mental data. [l3l In particular, the half-width at half- 
maximum at q = Q = 7r(l, 1) equals to Ag = 0.02l7r in 
Fig- [111(a) and Aq = O.OSStt in Fig.[T3](b) in comparison 
with the values Aq w 0.27r at T = 5 K and Aq w 0.247r 
at T = 70 K in Fig. 3 of Ref. ^ being about 7-8 times 
larger than our results. We believe that the difference 
in the damping may be explained by impurity scattering 
due to a disorder in the chains in the YBCOg.e sample or 
by an inhomogeneity produced by dynamical stripe fluc- 
tuations. The latter may be also the reason why we have 
not found a change of the dispersion topology above Tc 
observed in Ref. [l7[ which was related to the appearance 
of an electronic liquid-crystal state. 



IV. CONCLUSION 

A detailed study of the DSS in the superconduct- 
ing state has revealed the important role of the spin- 
excitation damping in the RM phenomenon. We have 
found that the low-temperature damping essentially de- 
pends on the gap wq ~ in the spin-excitation spec- 
trum at the AF wave vector Q, while the opening of a 
superconducting gap 2A(T) below Tc is less important. 
Since the energy of the RM E^- 5 fee 7c docs not show 
temperature dependence, at T < Tc, the spin gap at E^ 
plays the dominant role in the suppression of damping, 
since the superconducting energy 2A(T < Tc) <^ E^. 
This follows from Eq. (j25p for the self-energy, where in 
the decay of a spin excitation, besides a particle-hole 
pair, the contribution from an additional spin excita- 
tion is taken into account. This is in contrast to the 
particle-hole bubble approximation in Eq. (j26p which is 
used in the RPA (see, e.g., Refs. [13, [l!]) and in similar 
approximations in memory-function theories. [38l.l40ll42j 
In those approximations, the spin-excitation damping is 
much larger in the normal state and reveals a spin gap 
only at T <C Tc, where the RM appears (see Fig. [5]). The 
gapped spin-excitation spectrum in the full self-energy 
p5t greatly suppresses the damping at T = which re- 
sults in a comparable damping both in the superconduct- 
ing state (either of d-wave or s-wave symmetry) and in 
the normal state, as demonstrated in Figs. [T] and |31 The 
damping is strongly decreasing in the underdoped region 
(see Figs. [5] and [3]) bringing about a much stronger RM 
seen also above Tc (Fig. [T]). 

The weak temperature dependence of the RM in 
cuprates can be contrasted with the magnetic RM in 
iron-arsenide superconductors. In particular, a perfect 
scaling of the RM energy with the superconducting gap 
below Tc = 25 K observed in the BaFei.85Coo.i5As2 crys- 
tal [G^I supports the spin-1 exciton scenario expected in 
AF metals where the SDW instability determines the 
spin-excitation spectrum. Contrary to iron-arsenides, 
strong electron correlations in cuprates lead to a large AF 
exchange interaction between localized copper spins that 
results in temperature independent spin- wave-like excita- 
tions. A scenario of preformed Cooper pairs in cuprates, 
which may also explain a weak RM temperature depen- 
dence, seems to contradict experiments (for discussions 
see Ref. 

With decreasing hole concentration the energy of the 
RM decreases and shows the dependence close to that 
observed in neutron-scattering experiments, E^ ^ bk^Tc 
(see Fig. [11]). Due to the important role of gapped spin 
excitations in the damping of the RM, its energy E,- does 
not critically depend on the superconducting gap energy 
2A(T) and hence, on temperature and peculiarities of the 
electronic spectrum in cuprates, contrar y to the theories 
based on the RPA (see, e.g., Refs. 0, H^)- In partic- 
ular, Er is found to be larger than 2Ao as observed in 
experiments (see Figs . [T^ and [T5)) . 

It should be stressed, that the damping of the RM 
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in our microscopic theory within the t-J model is de- 
termined by the kinematical interaction induced by the 
kinetic energy t of electrons moving in a singly occupied 
Hubbard subband. This interaction is absent in con- 
ventional fermion models in which strong electron cor- 
relations are neglected and the spin-electron scattering 
is determined by a phenomenological interaction with a 
coupling constant as a fit parameter. 

Finally, let us note that our approach, using the MFA 
for the electronic spectral functions, Eqs. ([22l) and ([23l) . 
in the computation of the self-energy ([25]) . has to be con- 
sidered as a first step towards a fully self-consistent the- 
ory of the DSS in the t-J model. However, we believe 
that the consideration of more accurate fermionic GFs in 
Eqs. ([T6l) and ([TTl) beyond the Hubbard-I approximation 
(see, e.g., Ref. [65]) should not change our main conclu- 
sions, since the electron dispersion, as shown in Fig. [2l 
do not play an essential role in our theory contrary to 
the spin-exciton scenario. 



Acknowledgments 

We thank V. Hinkov for fruitful discussions. Partial fi- 
nancial support by the Heisenberg-Landau Program of 
JINR is acknowledged. One of the authors (N.P.) is 
grateful to the MPIPKS, Dresden, for the hospitality dur- 
ing his stay at the Institute, where a part of the present 
work has been done. 



Appendix A: Mode-coupling approximation 

To calculate the self-energy (|4]) we use the MCA for the 
time-dependent multi-particle correlation function which 
appears in the spectral representation of the relaxation 
function 



S"(q,c.) = , 



1 



-[/(q, -c^)-/(q,c^)], 



2cj m(q) 

/(q,c.) = / o^^e-*(^q|5^q(<))p™p-, (Al) 



where 

- St = [[St. {Ht + Hj)] , [Ht + Hj)] EE ^ Ft, (A2) 



determines the force correlation functions of the force op- 
erators F" denoted by the index a — tt, tJ, Jt, J J. As 
discussed in Ref. [i^, at a sizeable doping (5 '--^ 0.1 con- 
sidered in this paper, only the term Ff^ can be taken 
into account, since all other terms give negligible contri- 
butions. For this term we have 



Following the reasoning of Ref. [45|, in Eq. (jA3p only 
products of operators on different sites are taken into 
account. After the Fourier transformation to the q-space 
we obtain the force-force correlation function 



■[Fl't\F^\t)) = i2tr 



E E 

qi,q2 q'i,q^ 

K' O' n' H+ , ,t\ (A4) 
q3'i2^ii qi^qi^qs^ ' ^ ' 

Aq3,q2,qi ^qi,q2:q3 



'^qi:q2,q3 qi,q2,q3\ 

where qs = q — qi — q2 and q'3 = q — q'l — q'2- Here 
we introduce the vertex function 

Aqiq2q3 = '^{lm+<i2 ~ 7qi ) 793 + 7q2 - 7qi+q3 : (^5) 

where the terms linear in jq reflect the exclusion of terms 
in Ft with coinciding sites. 

In the MCA we assume that the propagation of 
electronic- and bosonic-type excitations at different lat- 
tice sites in Eq. (|A3P occurs independently which results 
in the decoupling of the correlation function (IA4[) into the 
corresponding single-particle time-dependent correlation 
functions. As it turned out by numerical evaluations (see 
also Ref. |45,]), the contribution from the charge excita- 
tions given by {X°° + XJ'^) in Eq. (|X3l) can be neglected 
in comparison with the spin-excitation contribution given 
by Xj — St ■ In this approximation we obtain the spin- 
diagonal correlation functions for the normal state 

(Ki,qi.qi]t|^f^„q.,q3W) 



= (^qt^Iq^^qi I K^A^K^ 

= (x'j^x^^it)) {sz^XM iKsKm 



X <^qi,q'j '^q2,q2'^q3,q3 



(A6) 

In the superconducting state we additionally take into 
account pair correlation functions which appear in the 
spin off-diagonal terms 



{[Ki 



j'l^q.,q2.q3W) 



qi.qbisJ i""qi.q2,q3v 



= -{x%^x°qi{t)) {sZqXM {x%,x^'M 



X '^qi,-q3 ^q2,q2'^q3,-q'i 



(A7) 

By substituting the MCA correlation functions (IA6p and 
(|A7[) in Eq. (|Aip we obtain the mult i- particle correlation 
function 



/CO 1 
-oo -^^ „^ 



■UK 



J — oo 



qi>q2 

Aq3,q2,qJ(^q:^q°(0)(^q3°^^:(i)) 
~Aqi,q2,q3Aq3,q2,qi (^-qi ^qi (^)) (^-qa "^qs (^))] ' 

CT 

Using the spectral representation for the time-dependent 
correlation functions (see, e.g., Ref. [s^) 



V(^)h(l/7r)]Im((A|i?})„, 



(A9) 
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where /(w) is the Fermi (Bose) function n{uj) ( N{uj)), after integration over time t we derive Eq. ^TE\\ 
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